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Abstract 

We present a general theory for the equilibrium structure of cylindrical 

tubules and helical ribbons of chiral lipid membranes. This theory is based 

on a continuum elastic free energy that permits variations in the direction 

of molecular tilt and in the curvature of the membrane. The theory shows 

that the formation of tubules and helical ribbons is driven by the chirality 

of the membrane. Tubules have a first-order transition from a uniform state 

to a helically modulated state, with periodic stripes in the tilt direction and 

ripples in the curvature. Helical ribbons can be stable structures, or they can 

be unstable intermediate states in the formation of tubules. 
PACS numbers: 87.10.-fe, 61.30.-v, 68.15.+e, 87.22.Bt 
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I. INTRODUCTION 



Chiral amphiphilic molecules can self-assemble into microstructures with a variety of 
morphologies. Two of the most interesting morphologies, both for basic research and for 
technological applications, are cylindrical tubules and helical ribbons Tubules are 

bilayer or multilayer membranes of amphiphilic molecules wrapped in a cylinder, as shown 
in Fig. [I] M. They are observed in a variety of systems, including diacetylenic lipids [|T],0, 
bile |J, surfactants M, and glutamates M. In diacetylenic lipid systems, the tubule diameter 
is typically 0.5 /im and the tubule length is typically 50-200 /zm; in other systems, the tubule 
dimensions are typically several times larger ||. In most of these systems, tubules exhibit 
a helical "barber-pole" pattern on the surface of the cylinder. Helical ribbons are similar 
microstructures, consisting of long twisted strips of membranes with their edges exposed to 
the solvent, as shown in Fig. |2| f9|. In some cases, helical ribbons are unstable precursors 
to the formation of tubules; in other cases, helical ribbons appear to be stable. Cylindrical 
tubules have been studied extensively for use in several technological applications, such as 
electroactive composites and controlled- release systems @. Helical ribbons have not been 
used in technological applications, but they have also been studied extensively as part of an 
effort to rationally control the self-assembly of tubules. 

There have been three general approaches to the theory of tubules and helical ribbons. 
First, de Gennes argued that a membrane of chiral molecules in any tilted phase must 
develop a spontaneous electrostatic polarization |10 . This polarization can induce a narrow 
strip of membrane to buckle into a cylinder. The original theory of de Gennes described 
buckling along the long axis of a cylinder, but a straightforward modification describes 
helical winding due to electrostatic interactions. This theoretical approach predicts that 
adding of electrolytes to the solvent should increase the radius of the resulting tubules, 
because the electrostatic interaction would be screened by electrolytes in solution. However, 
experiments have shown that electrolytes in solution do not affect the formation or radius of 



tubules [fTI], except for the particular case of tubules composed of amphiphiles with charged 
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head groups [T^|. Thus, electrostatic interaction is very probably not a dominant factor in 



tubule formation. 

As an alternative theoretical approach, Lubensky and Prost derived a general phase 
diagram for membranes with in-plane orientational order, which predicts cylinders as well 
as spheres, flat disks, and tori [ I3fl . Within the cylindrical phase, the cylinder radius and 



length are determined by a competition between the curvature energy and the edge energy. 
One important prediction of this theory is the scaling r oc i 1 ^ 2 between the tubule radius 
r and the tubule length t. However, experiments have found no correlation between r and 
I. Rather, in a typical sample of tubules, r is quite monodisperse while I varies widely |14 



Thus, the competition between curvature energy and edge energy is also not a dominant 
factor in tubule formation. 

A third approach, which is consistent with the experimental results on lipid tubules, is 
based on the chiral packing of the molecules in a membrane. Helfrich and Prost have shown 
that a chiral membrane in a tilted phase will form a cylinder because of an intrinsic bending 
force due to chirality JHJ. This bending force arises because long chiral molecules do not 
pack parallel to their neighbors, but rather at a nonzero twist angle with respect to their 
neighbors. If the molecules lie in bilayers, and are tilted with respect to the local layer 
normal, the favored twist from neighbor to neighbor leads the whole membrane to twist into 
a cylinder. 

Several investigators have generalized the original Helfrich-Prost concept of an intrinsic 
chiral bending force in different ways. Ou-Yang and Liu have developed a version of this 
theory based on an analogy with cholesteric liquid crystals [Tj|. Chappell and Yager have 



developed an analogous theory in which the direction of one-dimensional chains of molecules, 
rather than the direction of molecular tilt, defines a vector order parameter within the 



membrane ||17||. Nelson and Powers have used the renormalization group to calculate the 



effects of thermal fluctuations on tubules ||18|| . This calculation predicts an anomalous scaling 
of the tubule radius as a function of the strength of the chiral interaction. Chung et al. have 
considered the full elastic anisotropy of a membrane, and have related the pitch angle of 



tubules to a ratio of membrane elastic constants |J. They have also made predictions for 
the kinetic evolution of helical ribbons into tubules. 

In our earlier paper, we generalized the Helfrich-Prost concept in another way [[HJ. In 
that theory, we began with a three-dimensional (3D) liquid-crystal free energy and applied 
it to the case of chiral molecules in a membrane. We showed that the chiral term in this 
free energy has two simultaneous effects: it leads the membrane to twist into a cylinder, and 
it induces a variation in the direction of the molecular tilt on the cylinder. Through this 
calculation, we predicted that tubules can form periodic modulated structures characterized 
by helical stripes in the tilt direction winding around the cylinders, as shown in Fig. |3|. These 
stripes are analogous to the stripes seen in thin films of chiral smectic liquid crystals p0|-|22|, 
but in a cylindrical rather than a planar geometry. We argued that these stripes correspond 
to the helical substructure that is often observed on tubules. 

In this paper, we further extend this theoretical approach to provide a more unified 
and systematic model of tubules and helical ribbons. First, we combine our earlier theory 
with the Chung theory to obtain a theory of tubules with both elastic anisotropy and tilt 
modulation. We derive explicit expressions for the average direction of the tilt and for the 
direction of the modulation, in terms of the elastic coefficients, and we show that these 
directions are different. Next, we show that the same theory can also be applied to helical 
ribbons. This theory shows that helical ribbons can be stable microstructures for a certain 
range of parameters; they are not necessarily intermediate states in the formation of tubules. 
Finally, we investigate the effects of the striped tilt modulation on the detailed shape of 
tubules. Based on recent models of lipid bilayers p3| , p4| , we expect a modulation in the tilt 
to induce a modulation in the curvature. We predict the curvature modulation illustrated 
in Fig. |], which resembles the Pp rippled phases of lyotropic liquid crystals p5|j26[] . 



The plan of this paper is as follows. In Sec. II, we construct a free energy for tubules 
and helical ribbons, which includes elastic anisotropy and includes the possibility of tilt 
modulation. In Sec. Ill, we use this free energy to predict the radius and tilt direction for 
tubules with uniform tilt. The results of that section are in exact agreement with the theory 



of Chung et al. In Sec. IV, we consider tubules with a modulation in the tilt. We propose an 
explicit ansatz for that modulation, and then minimize the free energy over the parameters 
in that ansatz. We find that there is a first-order transition from uniform to modulated 
tubules. In Sec. V, we apply the same theory to helical ribbons. The theory predicts the 
range of parameters in which helical ribbons are stable. In Sec. VI, we investigate ripples in 
the curvature of tubules and helical ribbons. We derive an explicit expression for the ripple 
shape that corresponds to our ansatz for the tilt modulation. Finally, in Sec. VII, we discuss 
these theoretical predictions and compare them with experimental results. 



II. FREE ENERGY 

In this section, we construct a free energy for tubules and helical ribbons. In this free 
energy, we suppose that the membrane is in a fluid phase, which may have hexatic bond- 
orientational order but does not have crystalline positional order. This assumption is based 
on the theoretical argument of Nelson and Peliti, who showed that a flexible membrane can- 
not have crystalline positional order in thermal equilibrium at nonzero temperature, because 



thermal fluctuations induce dislocations, which destroy this order on long length scales |27 
This assumption is also supported by two types of experimental evidence. First, Treanor and 
Pace found a distinct fluid character in nuclear magnetic resonance and electron spin reso- 



nance experiments on tubules [28]. Second, Brandow et al. found that tubule membranes 



can flow to seal up cuts in the membrane from an atomic force microscope tip p9| . This 
flow indicates that the membrane has no shear modulus ||30|| . Further information on the 
membrane phase comes from Thomas et al., who found long positional correlation lengths 



of 0.068 to 0.135 jum in synchrotron x-ray scattering studies of tubules ||31|| . Taken together, 



all of these experimental results suggest that tubule membranes are in a highly correlated 
hexatic phase. 

In this paper, we do not attempt to derive the free energy from a 3D liquid- crystal free 
energy, as in our earlier paper [ITUJ. Rather, we use a 2D differential-geometry notation, 
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which is more general. In this notation, the curvature of a membrane is described by the 
curvature tensor K ab . The orientation of the molecular tilt is represented by the unit vector 
in. As shown in Fig. [5], fa is the projection of the molecular director n into the local tangent 
plane, normalized to unit magnitude. 

A membrane may consist of a single domain, in which the curvature K ab is low and the 
tilt direction fa varies smoothly. Alternatively, a membrane may consist of many internal 
domains separated by domain walls. Domain walls are sharp boundaries between domains. 
Across a narrow domain wall, fa jumps abruptly and K ab may become high. Here, we will 
consider the free energy of domains and domain walls separately. 

First, we consider a domain of low curvature and smoothly varying tilt. Within a domain 
of a membrane, the free energy can be written as 

1 -k {K a a f + l -K>m a m h K ac K c b 
+X UP m a m c K b a ^/ge bc 
-\ h sV9t ab D a m b - 1 K a a D b m b 



dA 



1 



a b 



(D a m° 



(2.1) 



Here, the surface area element is dA = y/gd 2 x, e ab is the antisymmetric symbol, and D a 
is the covariant derivative within the membrane. Throughout this work, we shall imply 
summation over repeated indices, such as a in the first term above. The first term of F is 



the standard Helfrich bending energy of the membrane |32| . The coefficient k is the isotropic 
rigidity. The second term represents the anisotropy of the rigidity. The coefficient k' is the 
difference between the energy required to bend the membrane parallel to the tilt fa and 
the energy required to bend it perpendicular to fa. In general, k' can be either positive or 
negative. The third term, introduced by Helfrich and Prost ||15|| , is a chiral term that favors 
curvature in a direction 45° from fa. The coefficient Ahp is a measure of the magnitude of 
the chiral interaction between the molecules in the membrane. The fourth term is another 
chiral term, which favors a bend in fa. This term was introduced by Langer and Sethna in 
the context of flat films of chiral liquid crystals |[20| . The coefficient Als should have the 



same magnitude as Ahp; both coefficients are comparable to Kq in the notation of our earlier 



paper |19j . However, there is no symmetry reason why Als and Ahp must be equal. The fifth 
term, with coefficient 7, is the coupling between curvature of the membrane and splay of 
in. It can be understood in two equivalent ways: (a) the curvature of the membrane breaks 
the symmetry between the two monolayers of a bilayer and hence induces a splay in m, or 
(b) a splay in m induces a curvature of the membrane to keep the molecules more nearly 
parallel in three dimensions. This term is equivalent to terms that have been considered in 
the Lubensky-MacKintosh theory of Pp> rippled phases |]23| , ^| . The final pair of terms is the 



2D Frank free energy for variations in m. It provides an energy penalty for tilt variations, 
and hence limits those variations. The coefficient c is a single Frank constant. 

We now consider a domain wall between two domains. Across a domain wall, the tilt 
direction m changes abruptly. The curvature K a b may also become large, giving a "crease" 
in the membrane. For those reasons, different parts of the molecules come into contact at 
a domain walls than inside the domains. A domain wall therefore costs an energy ew per 
unit length. The domain-wall energy cannot be derived from the free energy ( |2.1[ ) within 
a domain; rather, it is an additional parameter that describes narrow walls where that free 
energy does not apply. For any morphology that includes domain walls, with a total wall 
length Lw, the total free energy is F tot = F + ey/Lw- 

In the free energy (|2 . 1| ) , the term \\ J c,^J~ge a i ) D a m b is a total derivative. If the curvature is 
constant, as in a perfect cylinder, the term ^K^D^vr^ is also a total derivative. One might 
think that these terms integrate to constants and hence do not affect the morphology of a 
membrane. However, these terms are significant if the membrane consists of internal domains 
separated by domain walls, because they then become line integrals along the domain walls. 
Indeed, these terms favor the formation of a series of domains separated by domain walls, 
because they can contribute a negative free energy for each domain. Thus, in the rest of this 
paper, we will retain these terms, and we will explicitly compare their negative contribution 
to the free energy with the energy cost of the domain walls. 

The free energy ( p.l|) does not contain an exhaustive list of all the terms permitted 



by symmetry; other terms are also possible. For example, the term (m a m b K a b) 2 would 
give an additional elastic anisotropy, which has been considered by Chung et al. H. The 



combination {^fge a bD 



a m b\2 



(Db7n ) would give a difference in the Frank constants for 



bend and splay. Rather, this free energy is only intended to include the terms that generate 
distinct physical effects, while the omitted terms give higher-order details of the structure. 

If the membrane forms a perfect cylinder with radius r, the free energy ( p.l|) simplifies 
greatly. In the standard cylindrical coordinates {0,z), the curvature tensor becomes 

{ -l/r o\ 


The tilt director field can be written as rh(8,z) = (cos 0, sin 0). As shown in Fig. |5], (p(8,z) 
is the angle in the local tangent plane between the tilt direction and the equator of the 
cylinder. The free energy then becomes 



K 



\ 



(2.2) 
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— K ( — ) H — k' ( - ) cos 2 . 
2 \rJ 2 \rj 

— Ahp ( - ) sin cos 4> 



AlsV x m + 7 I -J V • m 



+ 2 C 



(2.3) 



where dA = rdOdz, 



V x rh = -de sin <p — d z cos < 

r 



(2.4) 



and 



V • m = -de cos0 + d z sin</>. (2.5) 

r 

Note that the 2D cross product defined by Eq. ( |2.4| ) is a scalar. This form of the free 
energy is generally more convenient to work with than the more general form ( |2.1| ). The 
interpretation of these terms is exactly as discussed above. In the following sections, we will 
use this free energy to investigate the structure of tubules and helical ribbons. 
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III. TUBULES WITH UNIFORM TILT 



As a first step in understanding the implications of the free energy ( j2.1|) , we consider 
tubules with a uniform tilt direction m. This model of uniform tubules is equivalent to the 
model of Chung et al. ||. We present it here in order to review their results in our notation. 

Suppose that a tubule of radius r has a uniform tilt fh = (cos0o,sin0 o ), at an angle 0o 
away from the equator, as shown in Fig. ^. The entire tubule is a single domain, with no 
domain walls. In terms of r and <fto, the free energy per unit area is 

F 1 f 1 ^ 1 / f 1 ^ 2 . 

— = -ft - H K - COS 0n 

-Ahp (~j sin 0o cos 0o- (3.1) 
Minimizing this expression over r and 0o simultaneously, we obtain 

«V4 (« + «/)V4 L.l/2 + ( K + K /)l/2" 



and 



(3.2) 



arctan 



'ft + ft'\ 1/4 ' 



(3.3) 



These results lead to two important conclusions. First, the tubule radius r scales inversely 
with the chirality parameter Ahp- In a nonchiral membrane, where Ahp — > 0, the radius 
r — > oo. The coefficient of 1/Ahp is a combination of the elastic constants k and k'. Second, 
the tilt direction 0o is determined by the ratio of the energy cost k + k' for bend parallel 
to the tilt direction over the energy cost k for bend perpendicular to the tilt direction. In 
an isotropic membrane, where k' = 0, we obtain O = 45°. In an anisotropic membrane, k' 
may be positive or negative, and hence 0o may be greater or less than 45°. 

Chung et al. originally derived Eq. ( ^.3[ ) in order to analyze experimental data on 
tubules in bile. In experiments on bile systems, they found two types of tubules. Both types 
of tubules show clear helical markings, but the direction of these markings are different: 53.7° 
from the equator in one type of tubule and 11.1° from the equator in the other. Chung et al. 



assumed that these helical markings are aligned with the molecular tilt, so that 0o — 53.7° 
and 11.1° in the two types of tubules. They then used Eq. ( |3.3| ) to determine the values of 
(k + k')/k that correspond to these values of <po. The resulting values are (k+k')/k = 3.4 and 
0.0015, respectively. The ratio of 3.4 is quite reasonable, but the ratio of 0.0015 is surprising. 
It seems unlikely that a bend in the membrane perpendicular to the tilt direction would cost 
almost 1000 times more energy than a bend parallel to the tilt direction. For that reason, we 
must re-examine the assumption that helical markings are aligned with the molecular tilt. 
In the following section, we propose an alternative interpretation of the helical markings, 
which can explain this anomaly. 



IV. TUBULES WITH TILT MODULATION 

In this section, we relax the assumption that the tilt direction m is uniform everywhere 
on a tubule. Instead, we allow m to vary as a function of the cylindrical coordinates 9 and 
z, with m(6>, z) = (cos0(#, z), sin <p(9, z)). We show that the ground state of the tubule can 
have a periodic helical modulation of m. 

As a specific ansatz, we suppose that there are n distinct stripes in the tilt direction. In 
the terminology of Sec. II, each stripe is a single domain of the membrane. The stripes are 
separated by n distinct domain walls. The stripes and domain walls run around the cylinder 
helically, as shown in Fig. [| We further suppose that the tilt direction <ft varies linearly 
across each stripe. Across each domain wall, <fi changes sharply, almost discontinuously on 
the scale of the stripes. The direction of the stripes and the domain walls is defined by 
the angle uj with respect to the equator of the cylinder, illustrated in Fig. [7[ Let 6' be the 
coordinate running along the uj direction, 

9' = 9 + -ta,nuj. (4.1) 
r 

In terms of this coordinate, our ansatz for the tilt direction can be written as 

0(M = 0o + ^(Vmod^V (4.2) 
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where 6' mod 2n/n runs from —ir/n to n/n. This ansatz has the sawtooth form shown in 
Fig. H It oscillates about the average value 0o with the amplitude A</>. Thus, our ansatz 
has five variational parameters: the radius r, the average tilt direction <fio, the amplitude of 
the tilt variation A0, the stripe direction u, and the number of distinct stripes n. In terms 
of these parameters, the stripe width is 

L= ( — ^)rcosc<j. (4-3) 
V n J 

We must now express the free energy ( |2.3|) in terms of the five variational parameters. 
To simplify the algebra, we suppose that A0 < 1; i. e., there is only a small variation in the 
tilt direction. We will discuss corrections to this approximation at the end of this section. 
The curl and divergence of the tilt then become 

nA^cos 

V x m = (4.4) 

2nr cos uj 

and 



V ■ m = . (4.5) 

27rr cos u 



With these expressions, it is straightforward to work out the gradient terms in Eq. 

The terms ^n'r" 2 cos 2 <fi — Ahp?" -1 sin cos require special attention. This combination 
of terms favors a particular orientation O of the tilt direction with respect to the equator 
of the cylinder. Because it gives an energy penalty for variations in the tilt away from the 
optimum direction, it depends on A0 as well as on 0o and r. We expand this combination 
as a power series in (<f) — <p ), then average it over the cylinder (i. e., average it over the 
coordinate^'). The result can be written as \k!t~ 2 cos 2 4>o — Anpr" 1 sin O cos^o + |^(A0) 2 , 
where 

1/2 

(4.6) 



1 

v = - 
6 



K \ / Ai 



+ 



VHP 



2r 2 J \ r 

The term ^i/(A0) 2 gives the energy penalty for variations of (p away from </>o- 

The domain-wall energy also requires special attention. As noted in Sec. II, the con- 
tinuum free energy of Eqs. ( p.l|) and (|2.3|) does not describe the narrow regions inside the 
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domain walls. Rather, we must explicitly add the domain-wall energy ew per unit length. 
In principle, ew can depend on O , A0, and uj. In this paper, we will neglect that possible 
variation and treat ew as a constant. In an area A of the membrane, the total length of 
all the domain walls is Lw = A/L. The domain- wall energy per unit area of membrane 
therefore becomes 

A 27rrcosc<j 

Putting all the pieces together, the total free energy per unit area for our ansatz now 
becomes 

Ftot 1 /1\ 2 1 ,/l\ 2 2 

" -K - + -K - COS 0o 



A 2 \rj 2 Vr 

-Ahp ( - ) sin 0o cos O + -z/(A0) 2 



-A 



LS" 



rj 2 
nA0 cos (0o — w) 7 nA0 sin (0 O — uj) 
2ttv cos a; r 2nr cos c<j 



+ 1/ nA0 \ 2 + e w n 
2 \2irrcosuJ 2ixr cos uj 

This expression must be minimized over the five variational parameters r, 0o> A0, u, and 
n. 

To do this minimization, we will make an approximation: We will treat n as a continuous 
variable rather than as an integer. This should be a good approximation for n>l, although 
not for n ~ 1. We will discuss corrections to this approximation, as well as the previous 
one, at the end of this section. After making this approximation, it is convenient to change 
variables from n to L = (2-nr cosu)/n. We also change variables from u to 5 = 0o — uj. The 
variable L is the stripe width, as noted above, and 5 is the difference between the average tilt 
direction and the stripe direction. In terms of these variables, the free energy ( |4.8|) simplifies 
to 

~A = 2 K [r) + 2 K It) °° S 00 



-Ahp (^j sin 0o cos 0o + ^z/(A0) i 
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A(f> cos 5 7 A(f> sin S 




(4.9) 



This expression for the free energy must be minimized over the five variational parameters 
r, 0o, A(j), 5, and L. 

As a first step, we minimize the free energy over the angle 5. We obtain 



This expression is equivalent to the corresponding result in our earlier paper This result 
is important because the angle 5 gives the difference between the stripe direction uj and the 
average tilt direction 0q. Because these two directions differ by an angle 5, we can explain 
the anomaly mentioned at the end of Sec. III. Recall that the experiments of Chung et al. 
found one type of tubules with helical markings at an angle of 11.1° from the equator. If these 
helical markings indicate the average tilt direction </> , then the ratio of rigidities (k + k')/k 
must be 0.0015, an anomalously low value. However, if the helical markings indicate the 
stripe direction, then this result for (k + k')/k does not hold. That ratio could be closer 
to 1, and the angle <po could be closer to 45°, while the stripe direction is 11.1°. Thus, our 
interpretation is that those experiments were measuring the stripe direction, not the average 
tilt direction. 

One might think that 5 would be very small because tan 5 scales as 1/r and r is large. 
However, in Sec. Ill we showed that r scales as k/Ahp- For that reason, tan 5 scales as 
(7Ahp)/(Als k )- We have already argued that the chiral coefficients Ahp and Als should have 
roughly the same magnitude. Hence, we obtain 



tan 5 



7 



(4.10) 



rA 



LS 



tan 5 



7 



(4.11) 



This ratio can be large or small. 



For the next step in the calculation, we substitute Eq. ( (4.10[ ) for 5 back into Eq. ( [4.9| ) 



for the free energy, to obtain 
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L tot 

~A 



2 \r. 



-K I - J COS n 



-A !U > ( sin 0o cos 0o + ^(A</>) 2 



1 

+ 2 C \ L 



1/2 



A0 



(4.12) 



Minimizing this expression over the stripe width L gives 

c(A0) 2 



(4.13) 



[AL + (7/r) 2 ] 1/2 A0- ew 
if this denominator is positive, or L — > oo otherwise. This expression for L is the same result 
that one would expect from studies of stripes in flat films |20^p2|. It is equivalent to the 



result for stripes in tubules in our earlier paper |19| , except that the dependence on A0 was 
not considered there. This expression shows that the Als and 7/r terms in the free energy 
both tend to induce stripes with a narrow width. If the energy gain from these terms, added 
in quadrature, exceeds the energy cost ew of a domain wall, then stripes will form; otherwise 
they will not form. If stripes do form, their width depends on the competition between the 
Frank free energy proportional to c, which resists variation in 0, and the Als and 7/r terms. 
We now insert Eq. ( |4.13| ) for L back into Eq. ( |4.12| ) for the free energy. The result is 



^tot 1 fl\ 2 , 1 



A 



-K 

2 \r 



+ ~K 

2 



cos <po 



-Ahp ( - ] sin 0o cos 0o + ^(A0) z 



1 

Yc 



\2 1 / 7 
A LS + 



1/2 



(4.14) 



if the quantity in brackets in the last term is positive. If not, then the last term is zero. 
Because only the last two terms in the free energy depend on A0, we will call them 
-F stripe (A0)/A We must now minimize -F s t r i P e over A0. This minimization cannot be done 
analytically, but it can be done graphically. In Fig. ^ we plot 

-^stripe & 

function of A0 for 

a sequence of values of v. The parameter u, defined in Eq. ( f4.6p , gives the energy penalty 
for variations of away from 0q. For large the only minimum of -F s tripe is for A0 = 0. 
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That is reasonable, because a large value of v suppresses all variations in 0. As v decreases, 
-^stripe develops a local minimum for A</> 7^ 0. At a critical value z/ c , there is a first-order 
transition from A<p = to A<p = A<ft c . For v < u c , the minimum at A<p 7^ becomes 
the absolute minimum of the free energy. To calculate v c and A0 C , we set -F str i pe = and 
<9F stripe /<9(A0) = 0. The result is 



1 



16ce^ 



A 2 , s + 



7 



(4.15) 



and 



A0 C = ^ m . (4.16) 

[^LS + (l/ r ) 2 } 



The corresponding value of the stripe width is 

Ace w 



C Ai s + ( 7 /r)2' 
and hence the gradient in the tilt direction <fi is 

A0 c _[AL + (7/r) 2 ] V2 



(4.17) 



(4.18) 



L c 2c 

In summary, for v > v c the tubules are uniform, with no tilt modulation. In that case, 
we have A0 = 0, L —>■ 00, and -F stripc = 0. At v = u c , there is a first-order transition from 
the uniform state to a modulated state with A0 = A0 C and L = L c . At this point, -F stri p e 
is still zero. For v < u c , the tubules remain in the modulated state. As v decreases, A(f) 
and L both grow larger and F s t r i P e becomes negative. It is not possible to derive a general 
analytic expression for Acf) and L in the modulated state. However, a good estimate would 
be Acj) m A(j) c and L w L c . 

We must now minimize the free energy ( 4.14 ) over the radius r and the average tilt 



direction 0o, which are the remaining variational parameters. For v > i/ c , we have -F s tripe = 0, 
and hence the total free energy is 

Ftot 1 /1\ 2 1 ,/l\ 2 2 
— r~ — —k - H — K - cos 0n 
A 2 \rj 2 Vr/ ^ 

-Ahp ( - ) sin0 o cos0 o - (4.19) 
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This is exactly the free energy considered in Sec. III. Hence, the solution given in Eqs. ( |3.2|) 
and ( p.3| ) still applies to the uniform state of tubules. For v < u c , we have -F stripe < and 
hence 

Ftot 1 fl\ 2 , 1 ,(1\ 2 2 , 

; K \~ ) + o K ~ cos 0o 



A 2 W 2 Vr 

-A HP sin 0o cos 0o + ^p- (4.20) 

Although we do not have a general expression for F str i pe , we see that -F str i P e can be expanded 
in a power series in r~ 2 . The constant term in that series does not affect the minimization. 
The term proportional to r~ 2 just gives a negative renormalization of the rigidity k. The 
higher-order terms in r~ 2 should be neglected, because we have already neglected terms of 
that order in our original curvature expansion for the free energy. For that reason, if we just 
use the renormalized rigidity in place of the bare rigidity k, then the derivation of r and 
0o in Sec. Ill still follows. We obtain Eqs. ( |3.2|) and ( |3.3| ) for r and 0o with the renormalized 
in place of n. 

At this point, we have minimized the free energy ( (4.9| ) over the five parameters in our 
variational ansatz: r, O , A0, 5, and L. We have shown that there is a first-order transition 
from a uniform state to a modulated state characterized by stripes in the tilt direction. We 
must now discuss corrections to two approximations in this variational calculation. 

First, we assumed that there is only a small variation in the tilt direction across a stripe: 
A0 <C 1. In some systems, A0 might be locked at a larger angle. For example, in a mem- 
brane with hexatic bond-orientational order, the angle A0 = 60° would be favored, because 
that angle would allow the tilt direction to jump from one bond direction to another bond 
direction across each domain wall. Such locking has been investigated in the context of flat 
liquid-crystal films by Hinshaw et al. |^TJ . In such a system, the transition from the uniform 
state to the modulated state would become more strongly first-order than we calculated 
above: the larger discontinuity in A0 would imply a larger latent heat of transition. Thus, 
the quantitative predictions of this section would no longer apply. Nevertheless, our quali- 
tative predictions would still apply, because tubules would have a first-order transition from 
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a uniform state to a modulated state whether or not A0 <C 1. 

Second, we did the variational calculation with the approximation that the stripe width 
L is a continuous variable, or equivalently, that the number of stripes n takes continuous 
rather than just integer values. In fact, because of the periodic boundary conditions going 
around a cylinder, the parameter n must be an integer, and hence L = (2irr cosu)/n is 
restricted to a discrete set of values. For that reason, there will be a series of jumps, or first- 
order transitions, between different integer values of n. Within each constant-n "phase," 
there will be one special point where the continuous minimization gives exactly the right 
answer: n = integer. At those special points, all the results of the minimization (for 5, L, 
A0, etc.) will be correct. Between those special points, there will be deviations from the 
predictions of the continuous approximation. The actual values should therefore fluctuate 
about the values from the continuous approximation. Thus, this approximation should give 
a reasonable estimate of the actual results. 

V. HELICAL RIBBONS 

In this section, we apply the theory for the modulated state of tubules, developed in the 
previous section, to helical ribbons. We draw an explicit analogy between a single stripe 
of a modulated tubule and a single ribbon. We show that a helical ribbon can be either a 
stable state of a membrane or an unstable intermediate state in the formation of a tubule. 
A stable helical ribbon has a particular optimum width. If more molecules are added to the 
ribbon, it becomes longer, not wider. By contrast, an unstable helical ribbon grows wider 
until it forms a tubule. 

The basis of our calculation is the geometry shown in Fig. |H]. A narrow ribbon winds 
helically around the z axis. We consider a linear variation of the tilt direction up to the 
ribbon edge. The tilt variation is in the direction specified by the angle u> with respect to the 
equator, as is the ribbon itself. Again, we describe the tilt variation using the 9' coordinate, 

9' = 9 + -ta,nuo. (5.1) 

r 
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which runs in the to direction. The ribbon edges are at 9' = ±9*, where 9* is a variational 
parameter. The maximum possible value is 9* = ir, at which point the edges of the ribbon 
collide, and the ribbon forms a tubule. Between the edges, our ansatz for the tilt direction 
is 

= 00 + ^0'. (5.2) 

This ansatz is a single period of the sawtooth wave of Fig. || As in Sec. IV, the ansatz has 
five variational parameters: the radius r, the average tilt direction 0o, the amplitude of the 
tilt variation A0, the ribbon direction to, and the angle 9*, which is the half-width of the 
ribbon in the 9' coordinate. The actual width of the ribbon is 

L = 29*rcosto. (5.3) 

We must now express the free energy (|2.3| ) of a helical ribbon in terms of these five 
variational parameters. The free energy is expressed per unit area of membrane, where the 
area does not include the gaps between the edges of the ribbon. For most of the terms, the 
free energy per unit area of a tubule, derived in Sec. IV, still applies. However, instead of 
the domain- wall energy ew, we must now use the edge energy €e per unit length for each of 
the two edges of the ribbon. This edge energy gives the energy cost of having different parts 
of the molecules exposed to the solvent at the ribbon edges. Hence, the total free energy 
per unit area in our ansatz for a helical ribbon becomes 

Ftat _ i m a , i , n\ 2 2 , 



-Ahp f -J sin 0o cos 0o + ^{^4>Y 

A0 cos (0o — to) 7 A0sin (0 O — to) 

LS 29*rcosto r 29*rcosto 
1 / A0 \ 2 e E 

+ 2 C \29*rcostoJ + 29*rcosto' ' 

This expression must be minimized over the five variational parameters r, 0o> A0, to, and 
9*. 
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At this point, we note that all five variational parameters are continuous variables; none 
of them is an integer, as in the previous section. Thus, we can immediately change variables 
from 9* to the ribbon width L = 26*r cosu), and from to to the difference angle 5 = (f)Q — uj. 
In terms of these variables, the ribbon free energy ( |5.4j ) simplifies to 

^tot 1 /1\ 2 1 ,/l\ 2 2 

~K I - I +-« hi COS 0o 



A 2 \rj 2 \r, 

-Ahp ( - ) sin O cos O + -z/(A0) 2 



-A 



LS" 



. ^ u ~7-u ' n 

r J 2 
cos S 7 A0 sin 5 



L r L 



1 /A0\ 2 2 ei3 
+ 2 C (~J (5 ' 5) 

Equation ( |5.5| ) for the ribbon free energy is exactly the same as Eq. ( |4.9| ) for the tubule free 

energy, but with 2e# substituted in place of ew However, for ribbons, all of the parameters 

really are continuous variables; that is not just an approximation. Thus, the theory of 

ribbons is mathematically simpler than the theory of tubules. 

Because the free energy for ribbons is equivalent to the free energy for tubules, we can 

carry over our results from the theory of tubules in Sec. IV. The direction of a ribbon, 

relative to the average tilt direction, is given by the difference angle 

tan<5 = 7 ~ ~- (5-6) 
rXhs k 



The ribbon width L is given by 



c(A0)- 
[A^ + (7/r) 2 ] 1/2 A0-2 e£ 



L = „. , , VT-i i/2 . - I (5-7) 



if this denominator is positive, or L —>■ oo otherwise. However, note that the maximum 
possible value of L is L max = 27rrcoso;. If L > L max , the edges collide and the ribbon 
forms a tubule. Thus, there is a stable ribbon with L < L max only for a certain window of 
parameters. 

For the amplitude A0 of the tilt variation, the analysis of Sec. IV applies again. At a 
critical value of u, 
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64ce| 



a ls + 



(5.8) 



there is a first-order transition from a state with A0 = and L — > oo, i. e., a state with no 
stable ribbon, to a state with A0 = A0 C and L = L c , where 

4:€ E 



A0 C 



[A!s + (7/r) 2 ] 



1/2 



(5.9) 



and 



8ce ; 



(5.10) 



C A^ s + ( 7 /r)2- 

If L c < L max , then a stable ribbon forms at this transition. Finally, the minimization over 
r and 0o goes through exactly as in the case of the modulated state of tubules. The results 
are the same as in Sec. Ill, but with the renormalized rigidity k,r in place of the bare rigidity 



K. 



VI. RIPPLES 

In this section, we consider the possibility of non-uniformly curved states of a tubule. In 
particular, we consider equilibrium undulated, or rippled, textures on otherwise cylindrical 
tubules. In part, this is because of the similarity of the model described above in Eq. ( p . 1|) 
to recent models of modulated phases of lipid bilayers p3| , |24|| and liquid crystal films [|33 



and the suggestion in Ref. [24| that chiral stripe textures such as those of Sec. IV may 



exhibit a rippled shape. There have also been unpublished experimental reports of regular 
periodic undulations of the surface of tubules |35[] . Quite generally, however, we expect that 



an underlying chiral stripe texture, with a spatially varying tilt field as described in Sec. IV, 
will lead to a shape modulation of the tubule, as shown in Fig. [|. This is because of the 
coupling of the molecular orientation, or tilt field, to membrane shape |23| . 

Given an underlying stripe texture of the molecular tilt rh(6, z), we study the correspond- 
ing tubule shape that results from the coupling of molecular tilt to membrane shape. Among 
the possible explicit couplings of molecular tilt to membrane shape, the term / yK%Dt > m b is 
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general to all membranes with in-plane orientational order. It is allowed by symmetry for 
both chiral and achiral lipid bilayers. Furthermore, the coupling constant 7 is expected to 
be of the same order as the bending modulus k |36|]. The second coupling in increasing 
powers of m and K, \^p^fgebjn a m c K^ is permitted only for chiral systems. Here, we shall 
consider a somewhat simplified model, in which we let k' = in Eq. fl2.1|) . We calculate the 
modulated shapes of tubules with tilt modulation in the limit that the ripple amplitude is 
small compared with the tubule radius r, which is of order 1 fim. 

For a rippled surface, of course, the curvature tensor K a b in Eq. (|2.1|) is no longer constant. 
In Sec. IV we described the tilt field fh(9, z) by a tilt angle 4>(9, z). Here, we consider also a 
ripple characterized by a deviation h(8, z) of the membrane surface away from a background 
cylindrical geometry. More precisely, the membrane position in three dimensions is given by 

X(6, z) = ((r + h) cos 9, (r + h) sin 9, z) , (6.1) 

where r is the fixed (average) radius of the tubule. In other words, we consider a small 
undulation h of the local cylinder radius. The following analysis is valid for h <r. 

We find it convenient to use coordinates (s,z), where s = r6. In these coordinates, the 
metric tensor for the tubule surface is 

'(1 + h/rf + {d s hf O s hd z h ^ 

K d s hd z h 1 + (d z h) 2 J 

where t a = d a X form a basis for the local tangent plane to the membrane. Through second 
order in the, presumed small, height modulation h, the surface area measure is dA = ^/gdsdz, 
where 

h\ / 1 



9ab — t a -tb 



(6.2) 



V^ = Vl det ^i- { 1+ rJ l 1 + 2 (V/l) J' (6 ' 3) 
and 

(V/1) 2 = (d s h) 2 + (d z h) 2 . (6.4) 

The curvature tensor is K a b = N ■ d a dbX, where N = t$ x t x /\t$ X t z \ is the surface normal. 
Through order h 2 , the mean curvature 
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K a =g K ab = -- + V h + — — 



(6.5) 



where g ab = {g a b) , and 



V 2 /i = d 2 h + d 2 z h. 



(6.6) 



We shall find that the tilt modulation described in Sec. IV leads to a modulation of h with 
the same period. For the stripe texture considered above, this period is the stripe width L. 
The last two terms in Eq. Q6.5|) are of order h 2 / (rL 2 ), which is smaller than V 2 /i ~ h/L 2 
by approximately h/r. Thus, we shall retain only the first two terms in Eq. Q6.5|) . 
To order h, the divergence of fh is given by 



D a m a = ( 1 — — | d s cos + d z sin + — — sin < 



(6.7) 



The first terms in Eq. ( |6.7| ) are of order A0/L. From Eqs. ( |4.18[ ) and (|5.6|) , we expect that 
this is of order 

1 



A0 ^ Als _ _ 
L c r 



(6.8) 



near the transition to the stripe texture, since we also expect that n and c are of the same 
order [Q. The last term in Eq. ( |6.7| ) is thus smaller than the other terms by a factor of 
approximately h/L. Below, we shall ignore the last term in Eq. ( |6.7| ). We comment on the 
validity of this below. 

The leading-order terms in Eq. ( |2.ip that depend on the ripple shape h are given by 



F h 



^fgdsdz 
dsdz 



K - {K a a f + 7 {KD (D b m b ) 



2 \r 2 



-7 [V 2 h 



*(^-2^ + (y 2 h) 



r 



(d s cos + d z sin . 



(6.9) 



+Ahp (sin cos [p 2 s h — d 2 hj + (sin 2 — cos 2 0^ d s d z h 



For a fixed background cylinder radius r and a (predetermined) modulating tilt angle 
which following Sec. IV can be described in terms of a single variable 8' given in Eq. ( [STT 
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the Euler-Lagrange equation for the height field h is determined by variation of Eq. ( |6.9|) 
with respect to h. To leading order in h/r, the result is an ordinary differential equation, 

Kd 4 s ,h - 7 <9 S 3 , cos (0 - u) + ^d 2 s , sin 2 (0 - u) = 0, (6.10) 

that can be integrated to yield 

ndg/h — 'ydg' cos (0 — uj) H — 5£ s i n 2 (0 — c<j) = constant, (6-H) 

where s' = rcosuoO' = scosu + zsinuj extends from —L/2 to L/2. Here, we have also used 
the fact that h is a periodic function. This does not allow, for instance, a linear term in 
Eq. ( 16.111) . For the tilt field given above in Eq. (f4.2|) , 



d s i cos (0 — uj) = — — sin (0 — uS) . (6-12) 

Lj 

Within one period of the modulated texture, the tilt angle can be represented as 

71 

- uo = 6 + A0-— s', (6.13) 
ZttL 

where again —L/2 < s' < L/2. For small A0, the solutions to Eq. ( |6.12j ) can be written 

h ^ hl m. a ) + ^ h2 m.,A, (6 . 14) 



where a and (3 are constants, and the coefficients h\ and hi depend on 7, k, Ahp, 8, 
and L. 

At this point, an additional assumption concerning the domain wall is necessary. Here we 
consider two possibilities, both of which, however, yield the same functional form (Eq. [|6.14|| ) 
for the height ripples on the surface of tubules. As noted above, we have not taken account 
of possible dependencies of the domain wall energy e-w on parameters of our model such as 
or uj. If we continue to assume as before that the domain wall is of infinitesimal thickness, 
and costs an energy ew per unit length that is independent of model parameters, which now 
include the possibility of a finite slope discontinuity of the membrane at the domain wall, 
then the constant in Eq. Q6.ll ) can be determined by minimizing the integrated free energy 
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of Eq. ( |6.9| ) using Eq. ( |6.11|) . Because of the presence of the domain wall, we must include 
total derivatives such as the term K,(V 2 h) jr in Eq. (|6.9|). This is equivalent to solving 



Eq. ( |6.11| ) with the constant on the right hand side set equal to 1/r. The result is given by 



Eq. ( ETT3D with 



and 

h 2 = -(^co S 5+^cos2s). (6.16) 

\0K L OK J 

For the above to be a periodic function in the range from — L/2 to L/2, we must have 

The other constant a is arbitrary. We note that in the above, the contributions to both hi 
and h 2 from the chiral and achiral couplings are of the same order, provided that Ahp ~ Als 
and 7 ~ k, since A(p/L ~ 1/r. 

On the other hand, if we assume a narrow but finite domain wall region (2), in which 
the elastic constants k, 7, and Ahp have the same values as in the region (1) above, and in 
which the tilt angle <fi varies again linearly, but in reverse to its variation in region (1), then 
the solution in region (1) can still be expressed by Eq. ( |6.14|) , where 



/i! = -^^sin5, (6.18) 

K L 



fi2 is given by Eq. ( |6.16| ), and (3 is given by Eq. ( p,17| ). This is valid for domain walls of 



width V <C L. The general result for L' « L is somewhat more complicated, although the 
general form of Eq. (|6.14| ) is still valid. In particular, periodicity of h in region (1) alone 



is no longer valid, and hence, the coefficient (3 differs from the value above. (A somewhat 
more general solution was derived in Ref. [[24]] for flat membranes.) In general, however, 
the membrane slope is continuous across the boundaries between regions (1) and (2) under 
the conditions of equal elastic constants k, 7, and Ahp- We note, however, that the chiral 
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and achiral couplings no longer contribute to the ripple amplitude at the same order. The 
dominant contribution to the ripple amplitude is from the achiral term / ~fK^Di J m b , while the 
contribution to Eq. ( |6.14[ ) from the chiral coupling Xnpy/geb c m a, m c K b is smaller by a factor 



of order A<p. 

For a stripe width L, the amplitude of the height modulation h is of order L 2 /r, where 
we have used Eq. ( |6.8|) . Thus, 

h /L\ 2 1 

-~ - ~—, 6.19 

r \r J n 2 

where, as in Sec. IV, n is the number of stripes on the cylinder. We have assumed that this 
is small. So, our analysis above is valid for L C r. In other words, we have calculated the 
shape of tubules with stripe textures in the limit that the stripe texture has both amplitude 
and period small compared with the tubule radius. Note also that h is smaller than L by a 
factor of L/r ~ 1/n. 



In Fig. 11, we show representative ripple shapes for various values of hi and h 2 . In 



Fig. |TT|a, we show the shape for hi < and h 2 = 0. This is the dominant contribution 
for small A0. This term is symmetric under s' — > —s'. The correction to this, smaller 
by order A0, is asymmetric. In Figs. |TT| b and |TT| c, we show the corresponding shapes for 
h 2 A(f) = 0.3/ii and h 2 A(p = hi. 

As a final point, note that the analysis in this section, as well as Sec. IV, applies to 
any membrane in a cylindrical morphology, regardless of how it formed. In particular, if 
any membrane is adsorbed onto a pre-existing cylindrical substrate (perhaps a microscopic 
wire or fiber), it can form stripes in the tilt direction and ripples in the curvature. These 
modulations can occur even if the membrane is not chiral: In a non-chiral membrane, they 
would be induced by the 7 term in the free energy, which couples membrane curvature to 
variations in the direction of the tilt. These modulations can reduce the free energy by 
concentrating the curvature into domain walls. 
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VII. DISCUSSION 



In this paper, we have presented a general theory of tubules and helical ribbons based 
on the concept of chiral molecular packing. This theory shows that tubules can have both 
uniform and modulated states. In the uniform state, tubules have a constant orientation of 
the molecular tilt with respect to the equator of the cylinder. In the modulated state, tubules 
have a periodic, helical modulation in the direction of the molecular tilt, and corresponding 
ripples in the curvature of the cylinders. In this section, we discuss the experimental evidence 
for these theoretical predictions. 

There are two types of experimental evidence supporting the concept that the formation 
of tubules and helical ribbons is due to chiral molecular packing. First, many experiments 
have seen helical markings that wind around tubules, giving tubules a chiral substructure. 
Clearly, helical ribbons always have a chiral structure. It is reasonable that the observed 
chirality of these microstructures results from a chiral packing of the molecules. Second, 
recent experiments have found that diacetylenic lipid tubules have a very strong circular 
dichroism, which indicates a local chiral packing of the molecules, regardless of whether 
a chiral pattern is visible on the surface of the cylinder |37|. The same diacetylenic lipid 
molecules in solution or in large spherical vesicles have very low circular dichroism. These 
results show that the molecular packing in tubules is chiral, while the molecular packing in 
spherical vesicles is not chiral. 

So far, there has not been any direct test of our prediction of tilt modulation — no ex- 
periments have been sensitive to the local direction of molecular tilt in tubules. However, 
the helical markings on tubules provide indirect evidence for this prediction. In some cases, 
these helical markings are boundaries between sections of tubules with different numbers of 
bilayers in the walls. However, in other cases, helical markings appear even when there is 
no detectable discontinuity in the number of bilayers, as in Fig. [l|. These helical markings 
are apparently stable, because they do not anneal away in time, and hence seem to be a 
characteristic of the equilibrium state of tubules. Our interpretation is that these helical 
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markings are the orientational domain walls predicted by our theory. In this interpretation, 
the domain walls are visible in electron micrographs because impurities accumulate there 
and colloidal particles from the solution preferentially adsorb there. Such preferential dif- 
fusion of impurities to orientational domain walls has been observed directly in Langmuir 
monolayers 



Of course, this interpretation of the observed helical markings provides only indirect 
support for our theory. For a more direct test of our theory, one would need an experiment 
that can directly probe the local direction of molecular tilt. One possible experimental 
technique is ffuourescence microscopy with polarized laser excitation. This technique has 
been used to observe variations in the local tilt direction in Langmuir monolayers [fjSfl . To 
apply this technique to tubules, one would either (a) use the intrinsic fluorescence of the 
constituent amphiphilic molecules, (b) attach a fluorescent group to the molecules, or (c) put 
a fluorescent probe into the membrane that forms tubules. One would then illuminate the 
tubules with polarized light from a laser source. Variations in the direction of molecular 
tilt would then lead to variations in the intensity of fluorescence, which could be detected 
using confocal microscopy or near-field scanning optical microscopy. Through this approach, 
optical techniques could detect the predicted helical modulation in the tilt direction in the 
modulated state. 

The ripples in tubule curvature predicted by our theory may be the modulations seen by 
Yager et al. |34| and by Thomas et al. . Electron micrographs taken in those experiments 



show very clear helical variations in the tubule curvature. However, other experiments have 
not observed any ripples on tubules, within the resolution of the electron micrographs. The 
preliminary data are not yet sufficient to explain why ripples are seen in some experiments 
but not in others. 

Our prediction of a first-order transition between the uniform and modulated states of 
tubules also has some experimental support. Nounesis et al. have measured the magnetic 
birefringence and specific heat of both single-bilayer and multi-bilayer tubules, as functions 
of temperature through the melting transition ||39|| . They find that single-bilayer tubules 
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undergo a second-order melting transition, with strong pretransitional effects, while multi- 
bilayer tubules undergo a first-order melting transition. Furthermore, single-bilayer tubules 
show an anomalous peak in the specific heat about 3 degrees below the main peak associated 
with melting into the untilted phase. This anomalous peak is consistent with our predicted 
transition between the modulated and the uniform states of tubules. The modulated state 
should occur in the 3 degree window between the anomalous peak and the main melting peak, 
where the membrane elastic constants are low, and the uniform state should occur at lower 
temperatures, below the anomalous peak. Thus, our theory may explain this heat-capacity 
anomaly. 

As a final point, we note that our theory for the modulated state of tubules leads to an 
interesting scenario for the kinetic evolution of flat membranes or large spherical vesicles 



into tubules. This scenario, first proposed in Ref. [^|, is illustrated in Fig. [12]. When a 
flat membrane or large spherical vesicle is cooled from an untilted into a tilted phase, it 
develops tilt order. Because of the molecular chirality, the tilt order forms a series of stripes 
separated by domain walls, as shown in Fig. |T2|a. Each stripe forms a ripple in the membrane 
curvature, and each domain wall forms a ridge in the membrane. Thus, the domain walls 
are narrow regions where different parts of the amphiphilic molecules come into contact with 
neighboring molecules and with the solvent. As a result, the domain walls become weak lines 
in the membrane, and the membrane tends to fall apart along those lines. The membrane 
thereby forms a series of narrow ribbons. These ribbons are free to twist in solution to 



form helices, as shown in Fig. 12b. Those helices may remain as stable helical ribbons, or 



alternatively they may grow wider to form tubules, as shown in Fig. |12| c. Note that this 
proposed mechanism for tubule formation can only operate if the initial vesicle size is larger 
than the favored ribbon width. If the initial vesicle is too small, it cannot transform into a 
tubule. This prediction is consistent with the experimental observation that large spherical 
vesicles (with diameter greater than 1 /im) form tubules upon cooling, while small spherical 



vesicles (diameter less than 0.05 /im) do not P0|PT[ ]. Thus, the theoretical prediction of 



stripes in the tilt direction gives some insight into the kinetics of the tubule-formation 
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process. 

In conclusion, in this paper we have shown the range of possible states that can occur 
in tubules. Tubules can have a uniform state, as was considered by earlier investigators, 
but they can also have a modulated state, with a periodic helical variation in direction of 
molecular tilt and in the curvature of the membrane. There is at least indirect evidence that 
the modulated state occurs in actual experimental systems. A more definitive test of this 
theoretical prediction requires direct experimental probes of variations in the molecular tilt 
direction. 
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FIGURES 




FIG. 1. Transmission electron micrograph of a tubule with adsorbed Pd/Ni catalyst particles 
on the surface ||]. The diameter of the tubule is approximately 0.5 /im. Reprinted from Ref. J19|| . 
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FIG. 2. Scanning electron micrograph of helical ribbons formed from a racemic mixture of D and 
L diacetylenic phospholipids. Note that the racemic mixture forms distinct right- and left-handed 
helical ribbons. The diameter of the helical ribbons is approximately 0.5 fim. Reprinted from 
Ref. p] with permission ©1988 Elsevier Scientific Publishers. 
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FIG. 3. Schematic view of the striped pattern in the tilt direction in the modulated state of 
a tubule. The arrows indicate the direction of the molecular tilt, projected into the local tangent 
plane. 
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FIG. 4. Schematic view of the ripples in the curvature in the modulated state of a tubule. The 
amplitude of the ripples has been exaggerated for clarity. 
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FIG. 5. Geometry of a tubule with radius r and length I. Here, n is the molecular director, 
m is the projection of n into the local tangent plane (normalized to unit magnitude), and (f> is the 
angle in the tangent plane between fa and the equator. 
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FIG. 6. Schematic view of a tubule with the uniform tilt direction m = (cos <fio, sin </>o), as 
indicated by the arrows. 




FIG. 7. Definition of the angle uj, which gives the direction of the stripes and the domain walls 
on a tubule with respect to the equator. 
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FIG. 8. Ansatz for the tilt direction (j)(9, z), in terms of the coordinate 9' = 6 + (z/r) tana;. 
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FIG. 9. Graphical minimization of the stripe free GI16rgy -^stripe ^ 

function of Acj), as discussed 

in the text. F s t r i P e is the sum of the two terms plotted in (a) and (b). (c) For large v, the only 
minimum is the uniform state with Acj) = 0. (d) As v decreases, another minimum develops for 
Acf> ^ 0. (e) At v = i/ c , there is a first-order transition from the uniform state (Acj) = 0) to the 
modulated state (Acj) = Acp c ). (f) For v < f c , the modulated state becomes the absolute minimum 
of the free energy. 
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FIG. 10. Schematic view of a helical ribbon, with a variation in the tilt direction across the 
width of the ribbon. The arrows represent the tilt direction on one side of the bilayer. 
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FIG. 11. Plots of the ripple shapes, as a function of the coordinate s' normal to the ripples. For 
small A(j), the dominant contribution is symmetric under s' — > — s', but corrections are asymmetric, 
(a) hi < and fr 2 = 0. (b) /i 2 A^ = 0.3/ti. (c) h 2 A(f) = h x . 
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FIG. 12. Scenario for the kinetic evolution of flat membranes into tubules, as discussed in the 
text, (a) When a membrane is cooled into a tilted phase, it develops stripes in the tilt direction, 
and then breaks up along the domain walls to form ribbons, (b) Each ribbon twists in solution 
to form a helix, (c) A helical ribbon may remain stable, or may grow wider to form a tubule. 
Reprinted from Ref. [|37| with permission ©1994 American Association for the Advancement of 
Science. 
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